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On  the  -use  of  singular  yield  conditions  and 
associated  flow  rules'^ 

By  William  Prager^ 

Abstract.  It  is  well  known  that  the  use  of  Trescats 
yield  condition  frequently  leads  to  a simpler  system  of  equ- 
ations, for  the  stresses  in  a plastic  solid  than  the  use  of 
the  yield  condition  of  Mises.  In  most  cases  where  Tresca’s 
yield  condition  has  been  used  for  this  reason,  the  flow  rule 
associated  with  the  Mises  condition  has  been  retained,  how- 
ever. Following  Koiter*  (1),  it  is  shown  that  further  sim- 
plification results  from  the  use  of  the  flow  rule  associated 
with  the  Tresca  condition.  The  reason  for  this  is  discussed 
in  connection  with  two  examples  concerning  the  finite  enlarge- 
ment of  a circular  hole  in  an  infinite  sheet  of  perfectly 
plastic  or  work -hardening  material.  The  second  example  is 
probably  the  first  non-trivial  case  in  which  a problem  of 
finite  plastic  deformation  of  a work-hardening  material  has 
been  treated  in  closed  form  by  the  use  of  incremental  stress- 
strain  relations. 


^ The  results  presented  in  this  paper  were  obtained  in  the 
course  of  research  sponsored  by  the  Office  of  Naval  Re- 
search (Contract  N7onr-35801) . 
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INTRODOCTl'ON 

The  mechanical  behavior  of  a perfectly  plastic  solid 
is  defined  by  the  yield  condition  and  the  flow  rule.  The  yield 
condition  specifies  the  states  of  stress  under  which  plastic 
flow  will  occur.  For  each  of  these  states  of  stress  the  flow 
rule  specifies  the  components  of  the  plastic  strain  rate  to 
within  an  arbitrary  common  factor.  This  positive  factor  is  re- 
quired by  the  assumption  that  the  solid  is  inviscid:  the  state 

of  stress  necessary  to  cause  a certain  type  of  plastic  flow  does 
not  depend  on  the  speed  of  deformation.  In  the  following,  a set 
of  strain  rate  components  that  are  defined  only  to  within  a 
common  positive  factor  will  be  said  to  specify  a "flow  mechanism". 

For  an  isotropic  perfectly  plastic  solid,  the  orien- 
tation of  the  principal  axes  of  stress  does  not  enter  Into  the 
yield  condition.  This  condition  therefore  assumes  the  form 

f(o1,d2,o3)  =o,  [ll 

where  the  function  f must  be  symmetric  in  the  principal  stresses 
dl?  d2>  d3*  ***  so^d  retains  its  Isotropy  during  plastic 

flow,  the  principal  axes  of  the  strain  rate  must  coincide  with 
those  of  the  stress.  The  flow  rule  therefore  reduces  to  a 
relation  between  the  principal  stresses  d2,  d^  and  the 
principal  strain  rates  e-j_,  e2 , e^. 

It  is  convenient  to  represent  the  yield  condition  [l] 
geometrically  by  considering  d^,  d2 , and  d^  as  the  rectangular 
coordinates  of  a point  on  the  "yield  surface".  The  flow  rule 
may  then  be  taken  to  specify  a direction  at  each  point  of  the 
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yield  surface,  the  three  direction  cosines  being  proportional 
to  e^,  e3« 

Mises  (2)  proposed  the  flow  rule  that  associates  with 
each  point  of  the  yield  surface  the  direction  of  the  exterior 
normal  of  this  surface  at  this  point: 


p •«  *e  - 9f 

ere2,e3  " W. 


[2] 


In  Eq.  [2]  It  has  been  assumed  that  the  sign  of  the  yield  func- 
tion f has  been  so  chosen  that  the  exterior  normal  of  the  yield 
surface  indicates  the  direction  of  increasing  values  of  f. 

The  idea  of  using  the  yield  function  f as  the  "plastic 
potential1'  proved  useful  for  the  formulation  of  flow  rules  for 
crystals  and  other  anisotropic  solids.  In  recent  years,  the 
concept  of  the  plastic  potential  has  gained  added  importance 
because  the  theory  of  limit  analysis  is  based  on  this  relation 
between  the  yield  condition  and  the  flow  rule  (see,  for  instance, 

(3),  (4)). 


The  flow  rule  [2]  presupposes  that  there  Is  a uniquely 
determined  exterior  normal  at  each  point  of  the  yield  surface. 
Not  all  yield  conditions  used  in  the  mathematical  theory  of 
plasticity  satisfy  this  regularity  requirement.  For  example, 
the  yield  surface  corresponding  to  Tre sea's  condition  of  con- 
stant maximum  shearing  stress  is  a regular  hexagonal  prism.  At 
a point  on  an  edge  of  this  prism,  the  exterior  normal  is  not 
defined,  and  there  arises  the  question  how  the  flow  rule  [2] 
should  be  modified  at  such  a singular  point. 

Consider  first  a point  P on  the  edge  formed  by  two 


adjacent  fiat  or  curved  faces  of  a yield  surface.  Each  of  these 
faces  has  a unique  exterior  normal  at  P,  and  each  of  these  nor- 
mals may  he  considered  as  representing  a flow  mechanism  that  is 
possible  under  the  state  of  stress  represented  by  P.  It  is  then 
natural  to  assume  that  other  possible  flow  mechanisms  can  be 
obtained  by  combining  these  two  fundamental  mechanisms  in  a 
linear  fashion  and  with  positive  coefficients. 

In  the  case  of  Tre sea's  yield  condition,  for  instance, 
the  faces  of  the  prismatic  yield  surface  lie  in  the  planes 

o1-d2=  + cJ,  = ± <32”cj3  = + <5,  [3] 

where  d is  the  yield  stress  in  simple  tension.  For  points  on 
one  of  the  adjacent  faces  d^  - d2  “ 0 and  d^  ~ °3  = ° 9 flov/ 
rule  [2]  furnishes 


•^:e2:e3  = 1:  - 1:0 

P+] 

^:e2;e3  - 1:0:  - 1, 

respectively.  The  first  flow  mechanism  represents  pure  shear 
in  the  d-p  d2  plane;  the  second  mechanism,  pure  shear  in  the 
d^,  d^  plane.  The  point  representing  the  state  of  simple  ten- 
sion d^  = d lies  on  the  common  edge  of  these  two  faces.  Accord- 
ing to  the  flow  rule  proposed  above,  the  flow  mechanisms  possible 
under  this  state  of  simple  tension  are  characterised  by 

ei:e2:e3  = 1:  - X:  - (1  - X) , [6] 

whore  0 <,  X < 1,  The  three  terms  on  the  right-hand  side  of  Eq, 
[6]  are  obtained  by  multiplying  the  corresponding  terms  of  Eqs, 
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[4]  and  [ 5l  by  X and  1 -X,  respectively,  and  adding, 

While  the  type  of  singular  point  just  discussed  is  the 
only  one  occurring  in  connection  with  Tresca’s  yield  condition, 
other  types,  such  as  a vertex  of  a polyhedron  or  a cone,  can  be 
treated  in  a similar  manner. 

At  first  glance  it  might  seem  that  this  modification 
of  the  flow  rule  [2]  at  singular  points  of  the  yield  surface 
would  greatly  complicate  the  mathematical  treatment  of  problems 
of  plastic  flow.  As  was  pointed  out  by  Koiter  (1),  however,  the 
contrary  is  trues  a considerable  simplification  results  from  the 
joint  use  of  Tresca's  yield  condition  and  the  associated  flow 
rule.  The  reason  for  this  will  be  explained  in  connection  with 
the  following  examples  which  will  also  demonstrate  that  the 
simplification  is  by  no  means  restricted  to  the  case  of  small 
plastic  deformations  considered  by  Koiter. 

TRESCA’S  YIELD  CONDITION  IN  PROBLEMS  OP  PLANE  STRESS 

In  problems  of  plane  stress  one  of  the  principal 
stresses,  say  <5 , vanishes,  and  the  yield  condition  can  be  rep- 
resented  by  a curve  or  polygon  in  the  c^  plane.  If  Tresca’s 
yield  condition  is  use!,  the  yield  polygon  is  formed  by  the  lines 
with  the  equations 

^1  “ d2  ~ ^ j i ^ , °2  ~ i ^ [7] 

obtained,  from  [3]  by  setting  = 0.  Figure  1 shows  this  hexagon 
which  is  .an  oblique  section  of  the  afore-mentioned  hexagonal 
prism.  The  axis  of  this  prism  passes  through  the  origin  0 
and  forms  equal  angles  with  the  positive  axes  of  0-^  and  0^ 
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(the  latter  being  normal  to  the  plane  of  Fig.  1). 

Consider  a generic  point  P on  one  side  of  the  yield 
hexagon,  for  instance  on  AB.  The  exterior  normal  of  the  three 
dimensional  yield  surface  at  P is  not  contained  in  the  plane  of 
the  figure.  The  projection  of  this  normal  on  the  plane  of  the 
figure  is  normal  to  AB,  however.  Thus,  the  ratio  e-j_:e2  ^or  ^9 
state  of  stress  represented  by  P may  be  obtained  as  the  ratio 
between  the  direction  cosines  of  the  normal  to  the  yield  hexagcsm 
at  P.  Once  the  ratio  e^:e2  is  known,  the  condition  of  incom- 
pressibility 

ex  + e2  + = 0 [8l 

yields  the  ratios  e^:e2:e2« 

All  states  of  stress  occurring  in  the  following  ex- 
amples will  turn  out  to  be  represented  by  points  on  the  side 
AB  of  the  yield  hexagon.  With  0 < \<  1,  the  modified  flow  rule 
then  furnishes  the  following  information. 

(a)  State  of  stress  represented  by  point  A: 

<31  = - d,  = d3  ~ 0 and  el:e2:e3  = C 9 i 

(b)  State  of  stress  represented  by  interior  point  of 
segment  AB: 

0-j_  ” °2  = * d3  = 0 and  e1:e2:e^  = -1:1:0.  [10] 

(c)  State  of  stress  represented  by  point  B: 

d2  = 61  ~ d3  = 0 and  e-L:e2:e^  = - \ : 1: -(1-X),  [lij 

If  the  yield  stress  a remains  constant  during  plastic 
flow,  the  solid  is  called  perfectly  plastic;  if  the  yield  stress 
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increases  during  plastic  flow,  the  solid  is  called  work-hardening. 
For  the  most  general  type  of  work-hardening  solid  the  form  of  the 
yield  condition  could  also  change  during  plastic  flow,  for  in- 
stance the  hexagon  of  Fig.  1 could  gradually  change  into  an 
ellipse.  This  would  correspond  to  a gradual  transition  from  the 
yield  condition  of  Tresca  to  that  of  Hises.  This  case  will  not 
be  considered  in  this  paper,  however.  It  will  be  assumed  that 
during  plastic  flow  the  yield  hexagon  remains  centered  at  the 
origin  and  merely  increases  in  slse.  In  the  following  the  yield 
stress  of  the  virgin  material  will  be  denoted  by  dQ  and  it  will 
be  assumed  that  the  rate  of  hardening  d is  proportional  to  the 
rate  D at  which  mechanical  energy  is  dissipated  during  plastic 

flow:  . 

d = aD  , [12] 

where  a is  a constant. 

For  the  three  cases  considered  above  D = d-^e^  + dgeg  + 
has  the  following  values: 


(a) 

D - - c5ei> 

Cl3l 

(b) 

D — ( d^  ~ do ) 6 — — oe^, 

IM 

(c) 

D — deg* 

[15] 

FIRST  EXAMPLE 

:,^|.he  first  example  consider  the  finite  enlargement 
of  a circular  hole  in  a sheet  of  a perfectly  plastic  material. 
This  problem  has  been  treated  by  Taylor  (5)  who  used  Tresca* s 
yield  condition  but  the  flow  rule  associated  with  Mises*  yield 
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condltion,  The  comparison  of  Taylor's  work  and  the  following 
analysis  reveals  the  great  simplification  achieved  by  the  use 
of  the  appropriate  flow  rule, 

A circular  hole  of  the  radius  a0  in  an  infinite  plate 
of  the  uniform  thickness  hQ  is  to  be  enlarged  to  the  radius 
a = 1*4  a0  application  of  a gradually  increasing  uniform  pressure 
to  the  edge  of  the  hole.  Since  finite  plastic  deformations  will 
be  considered,  elastic  deformations  will  be  neglected.  The 
material  at  a sufficiently  large  distance  from  the  hole  must  then 
be  treated  as  rigid  because  it  will  not  reach  the  yield  limit. 
Under  these  circumstances,  radial  displacement  of  the  material 
near  the  hole  is  made  possible  only  by  a thickening  of  the  sheet. 

In  the  elastic  part  of  the  sheet  the  radial  stress  or 
and  the  hoop  stress  dQ  are  inversely  proportional  to  the  square 
of  the  radius  r,  and  oQ  = - cp  > 0*  Thus,  the  state  of  stress 
at  the  elastic-plastic  interface  is  represented  by  the  center  of 
the  segment  AB  in  Fig,  1.  Just  inside  this  Interface  the  yield 
condition  and  flow  rule  are  therefore  given  by  Eqs,  [10]  with 
the  subscripts  1,  2,  3 standing  for  r,  9,  z,  respectively.  Since 
ez  = there  is  no  thickening  of  the  sheet  and  radial  displace- 
ment is  therefore  impossible.  The  material  just  inside  the 
elastic-plastic  interface  accordingly  remains  rigid  even  though 
the  stresses  dr  and  satisfy  the  yield  condition 

d - d ~ ~ o , fih] 

r 9 0 

The  condition  of  radial  equilibrium  in  this  rigid  plastic  oone  is 

3d  d-d 

r . r 9 _ 

W + 


r 


0. 


[17] 
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Prom  [16],  [17]  and  the  boundary  condition  °r  = - = " 6q /2 

at  the  elastic  plastic  interface  (r  = p0)  it  follows  that 


dr  = " do(2  + log 
Cq  = Cq(^  - log 


[18] 


in  the  rigid  plastic  zone,  The  hoop  stress  is  therefore  decreas- 
ing as  one  progresses  from  the  elastic-plastic  interface  into  the 
rigid  plastic  zone.  It  follows  from  the  second  Ea.  [18  ] that 
the  hoop  stress  vanishes  for 

r = dQ  / \/e  =0*606pQ.  [191 

This  radius  which  will  occur  frequently  in  the  following  work 
will  be  denoted  by  p • The  state  of  stress  at  r = p is  repre- 
sented by  the  point  A in  Fig.  1. 

Just  inside  the  circle  r = p,  the  circumferential 
strain  rate  eQ  must  vanish  on  account  of  the  rigidity  of  the 
surrounding  material.  The  radial  strain  rate  ep  need  not  vanish, 
however.  The  flow  mechanism  is  therefore  represented  by  the 
normal  to  the  side  AF  of  the  yield  hexagon.- 

There  are  now  two  possibilities  regarding  the  variation 
of  stress  and  flow  mechanism  inside  the  circle  r =p,  If  the 
point  representing  the  state  of  stress  moves  from  A towards  F, 
the  flow  rule  requires  that  eg  = 0,  This  means  that  the  radial 
velocity  v must  remain  zero  inside  the  circle  r = p,  because 
©g  = v/r.  The  sheet  would  therefore  remain  rigid  even  inside  the 
circle  r = p , If,  on  the  other  hand,  the  point  representing  the 
state  of  stress  remains  at  A,  the  greater  degree  of  freedom  in 


All-86  -10- 

the  choice  of  the  flow  mechanism  at  this  singular  point  allows  a 
plastic  deformation  of  the  sheet  inside  the  circle  r - p • 

Since  the  sheet  thickness  h will  not  remain  constant 


inside  the  circle  r = p , the  equation  of  radial  equilibrium  has 


the  form 


[20] 


Since  = - o , = 0 when  the  state  of  stress  is  represented 

by  the  point  A,  this  equilibrium  condition  may  be  written  as 


follows* 

5“(rh)  = 0. 


21 


While  it  is  convenient  to  use  the  terms  "velocity”, 
"strain  rate" , and  "rate  of  dissipation  of  mechanical  energy", 
it  must  be  remembered  that  the  plastic  solids  considered  here 
are  Inviscid#  Consequently,  the  flow  procoss  is  independent  of 
the  time  scale,  and  any  variable  that  Increases  monotonically 
during  the  flow  process  may  bo  used  as  a measure  of  "time”#  In 
the  following,  the  radius  p will  be  used  in  this  manner# 


If  the  radial  velocity  is  denoted  by  v,  the  radial  and 
circumferential  strain  rates  are  er  = 0v/9r,  and  = v/r#  The 
strain  rate  In  the  direction  normal  to  the  sheet  is  (l/h)Dh/dp  , 
where  Dh/Dp  = 9h/6p  + v3h/8r  denotes  the  "material"  derivative 
of  h»  The  condition  of  Incompressibility  requires  that  the  sum 
of  the  three  strain  rates  vanish: 


ISh 

h dp 


+ 


v dlj 


h 3r 


+ 


jfe  + * = 0. 

or  r 


[22] 


Equa  Equations  [21],  [22],  and  the  boundary  conditions 


o* 


[2£ 


h = h 


v = 0 


at  r = p 
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define  the  functions  h = h(r,p),  v = v(r,p).  From  Eq.  C 21 3 and 
the  first  boundary  condition  [23]  it  follows  that 


h = hQ  £.  , 
0 r 


[24] 


With  this  expression  for  h,  Eq,  [22]  reduces  to  a differential 
equation  for  v.  The  solution  satisfying  the  second  boundary 
condition  [23]  is 

[25] 


v * 1 - X 

P 


(Since  the  radius  p is  used  as  a measure  of  "time",  the  "velocity" 
must  be  dimensionle ss.  ) From  Eq,  [25]  the  strain  rates  are  ob- 
tained as  follows 

[26] 


e = M = - 1 
er  9r  p» 


efl  = 2 = 1 . 1 . 
w r r p 


If  the  material  that  at  the  "time"  p is  at  the  radius 
r was  initially  at  the  radius  rQ,  conservation  of  mass  requires 
that 


2* 


PP 

Jr 


2 p 

hr  dr  = nh  (p  - rrt). 


[27] 


Substitution  of  [241  into  [27] yields 


r = 


2 2 
P + rX 


2p 


[28] 


Equation  [24]  can  therefore  be  written  in  the  form 

2 


h = 2h 


o ”2“  ~2 

P + ro 


[29] 


According  to  [9],  the  solution  just  obtained  will  be 
valid  only  if  e-j_:e2  lies  between  -1:0  and  -1:1,  The  first  of 
these  bounds  is  attained  at  r = p where  = 0 by  the  second  Eq, 
[263,  The  other  bound  is  attained  at  r = p/2  where  Eqs.  [26] 
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fur  nish  er  = 1/p,  = 1/p.  Equation  [28]  shows  that  r = p/2 

corresponds  to  r0  = 0.  Thus,  the  solution  remains  valid  even  in 
the  case  where  a pin  hole  is  enlarged  to  some  finite  radius. 

Formally  the  results  obtained  above  agree  with  those 
furnished  by  the  "simplified  analysis"  which  Taylor  (5)  attri- 
butes to  Bethe.  Taylor  rejects  this  analysis  because  the  ratios 
of  the  principal  strain  rates  vary  while  the  ratios  of  the  prin- 
cipal stresses  are  constant.  As  has  been  shown  above,  Bethe's 
results  can  be  obtained  from  a consistent  theory,  and  only  experi- 
ment can  decide  whether  the  behavior  of  a given  material  agrees 
better  with  the  formulas  of  Taylor  or  Bethe. 

SECOND  EXAMPLE 

The  second  example  differs  from  the  first  one  only  in 
so  far  as  the  sheet  material  is  now  supposed  to  be  work-hardening 
in  accordance  with  Eq.  [12],  Since  no  hardening  takes  place 
where  the  material  remains  rigid,  the  stresses  outside  the  cii  e 
r =p  are  the  same  as  in  the  previous  example.  To  keep  track  of 
the  progressive  hardening  of  the  material  inside  the  circle 
r ~ P » it  will  be  necessary  to  use  the  Lagrangian  coordinate  r0 
as  one  of  the  independent  variables,  the  other  independent  varia- 
ble being  p as  before. 

Again  cQ  = 0 inside  the  circle  r = p*  Since  or  = - d 
is  no  longer  constant,  however,  the  condition  of  equilibrium  [21] 
must  now  be  replaced  by 

£(rh0)  = T5r)ar0)  W-Jrhe)  = °- 


[30] 
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The  solution  of  this  equation  that  satisfies  the  boundary  con- 
ditions h = hQt  g - cQ  for  r = p is 

rho  = ph0<JQ  , [31] 


In  accordance  with  Eqs.  [12]  and  [13],  the  rate  of 
strain  hardening  is 


where 


[32] 


= = log  1 1 • [33] 

With  the  use  of  this  expression  for  e-p  Eq.  [32]  can  be  written 
as  follows: 


^-(<,1/a  1 1 ) = 0. 

8p  1 8rn  1 


The  circumferential  strain  rate  is 


[3*0 


e = v=l^  = llogr, 
6 r r ap  dp  & ’ 


[351 


and  the  strain  rate  in  the  direction  normal  to  the  plane  of  the 
sheet  is 


ez  = h * T-  los  h. 

z h 3p  Sp 


[36] 


The  condition  of  incompressibility  can  therefore  be  written  in 
the  form 

^Chr  | $£-  | ) =0.  [37] 

Those  solutions  of  the  differential  equations  [3*+]  and 
[37]  that  are  compatible  with  Eq,  [31]  and  satisfy  the  boundary 
condition  6 s cQ  for  r = p are 
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and 


r - 


. p a/ ( 1+cc ) 

tl  + a * «uV2+a)Al+a)] 

2 + a p 


[33] 


[39] 


As  before,  the  solution  represented  by  Bqs.  [31], 

[38],  and  [39],  is  valid  only  as  long  as  the  ratio  er:eg  remains 
below  -1.  Brora  Eqs*  [33]  and  [3?]  it  is  found  that  this  is  the 
case  as  long  as  r > p/2*  According  to  Eq.  [39]  this  means  that 
the  following  inecualities  must  be  fulfilled  if  this  solution 


is  to  be  valid: 


its 

> (£>2+a 


< 2-1-J  U + <1  + a)  $ KS) 


1+a 

2+a 


[40] 

[41] 


To  indicate  the  manner  in  which  Bethe’s  solution  for 
a perfectly  plastic  sheet  is  modified  by  work-hardening , assume 
that  a = 0.  5*  For  simple  tension,  this  corresponds  to  one  half 
of  a percent  increase  in  yield  stress  for  a strain  of  one  per 
cent.  Figure  2 shows  the  thickness  distribution  plotted  versus 
r/p  and  r/rQ.  The  solid  curve  corresponds  to  a =0.5  and  the 
dotted  curve  to  a =0,  It  is  seen  that  even  such  a small  amount 
of  strain  hardening  has  a significant  influence  on  the  thickness 
distribution*  According  to  Ec.  [4l]  the  solution  represented 


by  the  solid  curve  in  Fig.  2 is  valid  for  a/aQ  < 1.48. 


CONCLUSION 

It  has  been  shown  that  the  joint  use  of  Tre sea's  yield 
condition  and  the  associated  flow  rule  can  lead  to  a considerable 
siraolification  of  the  mathematical  work  in  problems  of  plane 
plastic  stress.  The  reason  for  this  is  the  following: 


A1  -86 


-15- 


corresponding  to  the  sides  and  the  vertices  of  the  yield  polygon, 
one  has  either  a unique  flow  mechanism  and  a one  parameter  family 
of  principal  stress  values,  or  unique  values  of  the  principal 
stresses  and  a one  parameter  family  of  flow  mechanisms*  In 
either  case  considerable  simplifications  result,  and  the  complete 
solution  is  built  up  of  zones  in  which  one  or  the  other  type  of 
simplification  pertains. 
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